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ABSTRACT 



We present a method to constrain cosmic topology from the distribution of astronomical objects projected on the celestial sphere. 
This is an extension of the 3D method introduced in Fujii & Yoshii (2011) that is to search for a pair of pairs of observed objects 
(quadruplet) linked by a holonomy, i.e., the method we present here is to search for a pair of celestial sphere n-tuplets for n > 3. 
We find, however, that this method is impractical to apply in realistic situations due to the small signal to noise ratio. We conclude 
therefore that it is unrealistic to constrain the topology of the Universe from the celestial distribution, and the 3D catalogs are necessary 
for the purpose. 
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1. Introduction 

Determining the shape of space is one of the most important 
themes in modern cosmology. Together with the assumption of 
the cosmological principle, Einstein's General Relativity con- 
nects a curvature of space with the average energy density of 
the universe, and recent observations suggest a flat A-CDM uni- 
verse (e.g., Otot = 1 .0050!°°JJ^'|' from WMAP-i-BAO-i-SN data, by 
Hinshaw et al. 2009). 

However, while General Relativity describes local geom- 
etry (curvature), it predicts nothing about global property of 
space, namely, topology of the universe. Global topology of 
the Universe is often assumed to be simply connected with no 
evidence, though it can be rather multiconnected. A multicon- 
nected space with a nontrivial topology can be imagined as a 27*^- 
polyhedron, called Dirichlet domain, whose K pairs of faces are 
glued mathematically by holonomies. An object passing through 
one face immediately returns through the glued face. As a result, 
multiple images of single objects, often referred to as "ghosts", 
appear like those in a kaleidoscope (for detail, see, e.g. Lachieze- 
Rey and Luminet 1995). 

Many methods for constraining cosmic topology are based 
on this prediction, e.g., circles-in-the-sky method (Cornish et al. 
1998) is to search for intersections of the last-scattering surface 
and the faces of our Dirichlet domain. They are circles with 
the same temperature fluctuation pattern in the CMB map, be- 
cause they are copies of a physically identical region. Various 
authors have searched for matched circles using the WMAP 
data, and obtained diverse results. For example, Aurich (2008) 
has found a hint of a 3 -torus topology, while Roukema et al. 
(2008) have found a signature of a Poincare dodecahedral space 
topology. Contrary to them, Cornish et al. (2004), Key et al. 
(2007), and Bielewicz & Banday (2011) have found no topo- 
logical signatures and obtained the lower limit of the size of our 
Universe. The most recent constraint is ~ 27.9 Gpc by Bielewicz 
& Banday (2011). These disagreements suggest an existence of 



methodological problems and bring our interests to the 3D meth- 
ods, i.e., those that are to use distributions of galaxies, galaxy 
clusters, or active galactic nuclei. 

3D methods are to search for unusual positional patterns of 
objects in a given astronomical catalog, due to an existence of 
topological copies (e.g., Lehoucq et al. 1996; Roukema 1996; 
Uzan et al. 1999; Marecki et al. 2005; Fujii & Yoshii 2011). 
For example, if we live in a 3-torus universe, a pair separation 
histogram (PSH) of a given catalog will show sharp spikes at the 
separations corresponding to its fundamental lengths (Lehoucq 
et al. 1996). Fujii & Yoshii (201 1) have introduced a new method 
that is much more sensitive to topological signatures than the 
prior ones. This method can apply to any of 17 multiconnected 
Euclidean spaces, even if the whole universe is comparable to 
the observed region in size, i.e., in a case that there are only a 
few topological ghosts. 

Our 3D method and others all require spectroscopical obser- 
vations to calculate the radial distances of objects. A wide field 
spectroscopic survey is usually difficult, especially for high red- 
shift objects in which we are most interested. Mainly in order to 
explore CMB data, Bernui & Villela (2006) have extended the 
PSH method to the pair angular separation histogram (PASH) 
method that does not require spectroscopic data when applied to 
astronomical objects. Their method, however, is so insensitive to 
topological signatures that the signal does not appear unless we 
average a number of different but statistically equivalent PASHs, 
which is almost impossible in practice. 

A significant number of high redshift objects (galaxies or 
quasars) are recently found, though most of them are not spec- 
troscopically observed yet. This is a good point in time to ex- 
amine whether these objects are usable for constraining cosmic 
topology or not. In section 2, we introduce another 2D method to 
constrain cosmic topology from the distribution of objects pro- 
jected on the celestial sphere, which is an extension of the 3D 
method of Fujii & Yoshii (201 1). In section 3, we show that the 
method is in principle applicable, but in a realistic situation it is 
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no longer useful, similarly to that of Bernui & Villela (2006). We 
conclude therefore that spectroscopically observed 3D catalogs 
are necessary, in order to reveal the topology of the Universe 
from the distribution of astronomical objects. Throughout the 
paper we consider flat universes with zero curvature. 

2. Method 

In this section we describe the methodology forjudging whether 
a given pair of n-tuplets (2n-tuplet) are linked by some holon- 
omy or not. Our assumption is that the universe has zero curva- 
ture (Euclidean geometry), and the distances from us to the 2n 
objects are unknown while their celestial positions are known. 

A convenient way of writing the holonomies is to use a 
4D coordinate system (w, x,y,z) where the simply connected 3- 
Euclidean space is represented as a hyperplane w = 1 (see, e.g., 
Fujii & Yoshii 201 1). Then every holonomy y in a flat universe 
can be written as a 4D matrix, y = yxTNT. where yx and 7nt are 
a translational part and a nontranslational part, respectively. If a 
pair of n-tuplets (JCi, • • • ,x„) and {x[, ■ ■ ■ ,x'J are linked by y, 
we have 

x\ = yxi = yTTNT-"^! = 7mXi + L, 

(1) 

where L is the translational vector. Note that a vector X is 4D, 
Z = (1,X, y,Z), but the 3D part {X, Y,Z) is important. We know 
all the mathematical possibilities for y^ (an identity, an n-th 
turn rotation for n = 2, 3,4, or 6, or a reflection), so these 3n 
equations have 2n + 3 unknown quantities: the distances from us 
to the 2n objects and L. If the distances are known from spectro- 
scopic observation, then the unknown quantities are just three, 
L = {L\,L2,Ly), so we derive 3n - 3 conditions as we did for 
n = 2 in Fujii & Yoshii (2011). In that work we introduced a 
new method that is to search for quadruplets satisfying the con- 
dition 

x\ -x\^ ywiixx -X2), (2) 

and showed that the method is extremely sensitive to topological 
signatures in a given catalog. 

In the case considered here, however, we do not know the 
distances and have to take another way. For this, we first denote 
the 3D positions of objects as follows: 
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where r is the unknown distance from us to the object, and the 
vector R - (i, y, z) represents the known celestial position such 
that x} +y^ +z!' = 1. For the simplest case where n = 3 and 
yNT = id, i.e., y is a translation, the following relations hold by 
eliminating L: 

nRi - r2«2 = r[R[ - r'^R'^, (3) 
nRi - rsRi = r[R[ - r'^R'^, (4) 
Then eliminating r2, ^"3, r'2, and r'^ from these equations gives 

[(Ri X R2\r2 + {Ri X R'^Xzi + iR'2 X Ri\z2]ri + 
[(/?; X R'^Ui + (R2 X R2)zz\ + (Ri X R'llz^Vi = 0, (5) 

[(Ri X R3\z'^ + (R3 X R'^lzA + {R'j X Ri\z3Vi + 
[(/?; X R'^),Z3 + iR'3 X R3Vz[ + (Rs X R[)zz'3]r[ = 0. (6) 



These are simultaneous equations for ri , r[ and should have a 
solution (ri, r[) (0, 0), so the following relation must hold: 



AD-BC = 0, (7) 

where 

A = (/?i X R2\Z'2 + {Rl X R'2)zZl + {R'2 X Rl)zZ2, (8) 

B = {R\ X R'^^h + {R'2 X R2)zz'i + {Rl X /?'i),z^, (9) 

C = (/?i X R3\z'3 + {R3 X R'3\Zx + {R'3 X Ri\Z3, (10) 

D = (/?', X R'^),Z3 + {R'3 X /?3),-z'i + {R3 X /?',),Z3. (11) 



This relation is written by the celestial positions of the objects, 
and thus spectroscopic observations are not needed to use it. 
Though we have considered here the special case that y is a 
translation. A, B, C, and D for other cases are calculated simi- 
larly. If y is a half-turn corkscrew motion or a glide reflection. 



A = -{Ri X R2)zZ2 + {Rl X R'2)zZi + {R'2 X Ri)zZ2, (12) 

B = {R\ X R'^Ui - {R'2 X R2W1 - {R2 X R'^\z'2, (13) 

C = -{Rl X R3\z'3 + {R3 X R'^)zZi + {R'3 X Ri\z3, (14) 

D = (/?; X R'^\z3 - {R'3 X R3)zz'i - {R3 X R'Ozz'3. (15) 



And if y is an n-th turn corkscrew motion for n = 4, 3, or 6, then 

A = (Rl X R2U2 + {{Rl X R2)zil + («2 X «l)zZ2)c + IR2 X («2 X ROI^i, (16) 

B = (R\ X R^).Z2 + {(«2 X ^i)zi'i + («2 X R\),z'2]c - \Ri X (Rj x (17) 
C = (ffi X ff3)zZ3 + {(R3 X R!^),zx + {R!^ X Ki).Z3)c + x (R3 x Ri)|,*, (18) 
D = (ff'i X R'jUi + l(«3 X R3),z[ + (R3 X R\kz'j]c - {R3 x (R^ x R[)ks, (19) 

where c = cos{2n/n) and s = sm{2n/n). Note that the cases for 
n = 1 and n = 2 corresponds to that of translation and half-turn 
corkscrew motion, respectively. 

Given a catalog with objects, we count the number of 
sextuplets (2n-tuplets for n = 3) satisfying the condition of 
AD - BC - 0, for each type of holonomies, but within a cho- 
sen tolerance s, i.e., \AD - BC\ < s. The number of such sextu- 
plets will be larger for a multiconnected space than stochastically 
expected for a simply connected one. We assign each object Xi 
an integer s,, the number of sextuplets that satisfy the condition 
\AD-BC\ < E and also include Xj as one of their members. An 
histogram for a multiconnected space will have some bumps in 
the large-.?, region, since topological copies contribute to signal 
more frequently. 

Similar calculations can be done for n > 4 where we have 
more additional conditions, however, the calculation time is 
roughly proportional to nC„ = N\/{N - «)!«!, and the calcu- 
lation for a large n is unrealistically time-consuming. In the next 
section we show the results for n = 3, the minimum value above 
which this method can apply, and discuss the possibility that we 
can constrain the topology of the Universe by the method. 

3. Simulations and discussions 

In order to see the applicability of the method described in the 
previous section, we generated toy catalogs in simply and multi- 
connected Euclidean spaces. As multiconnected spaces, we con- 
sidered the following six cases: 

- one with a pair of translations, 

- four with a pair of n-th turn corkscrew motions for n = 

2, 3, 4, and 6, 

- and one with a pair of glide reflections. 
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For each simulation, the obsei-ved region is a unit sphere cen- 
tered at the observer's position, and the translational distance is 

L= |L| = 1.4, which implies that l-£'"^;r(l-xV-^/T ~ 12% 
of the observed region is a copy of some part of the "first" copy 
of the universe. Now consider L to be c/Hq ^ 4.2 Gpc, where c 
and H() are the speed of light and the Hubble parameter, respec- 
tively. With this example scale, the radius of the observed region 
is ~ 3.0 Gpc, corresponding to z ~ 0.9. The effects of global 
inhomogeneity that had been investigated in detail in Fujii & 
Yoshii (2011) were not considered here: the axes for the n-th 
turn corkscrew motions and the reflectional plane for the glide 
reflections pass through the observer, and we chose the correct 
coordinate systems. 

Our toy catalogs were made as follows: we uniformly dis- 
tributed 50 real objects, and then their copies were generated 
for each multiconnected space. The total number of objects 
(50 real ones plus their copies) is different from simulation to 
simulation, since some of the copies were beyond the observed 
region. Next, another catalog in the simply connected space with 
the same number of objects was made for comparison. All ob- 
jects were projected on the celestial sphere, and sextuplets satis- 
fying the condition \AD - BC\ < e were searched for 

First we succeeded in detecting topological signal with the 
choice of e = 10 **. The results are given in Table [1] The number 
of sextuplets that satisfy the condition of \AD - BC\ < s, for each 
case, is given as "topological index". [Note that we had normal- 
ized the topological indices with respect to the total number of 
quadruplets in Fujii & Yoshii (201 1), but not here.] As expected, 
topological index for a multiconnected space is larger than that 
for the simply connected one, which successfully distinguishes 
the two topologies. The difference can more clearly be seen in 
s,-histograms as shown in Figure [T] Some bumps are seen in the 
large-s; regions for multiconnected cases, because of the topo- 
logical copies that contribute to the signal more frequently than 
a stochastic expectation. 

These results show that this 2D method in principle con- 
strains the topology of the Universe, however, the practical appli- 
cation of the method is not so straightforward since the condition 
considered here is optimistic and somewhat arbitrary: 

- The tolerance as small as e = 10 cannot be chosen when 
correctly considering the limited resolution, e.g., with the 
typical ground-based optical instruments, < 1 arcsec, which 
implies the uncertainties in s of the order of 10"^. Moreover, 
even if we have an ideal instrument, due to the peculiar ve- 
locity, the quantity \AD - BC\ will be larger than 10"^ even 
for a really topological sextuplet. For example, the typical 
peculiar velocity of v = 500 km/sec and the time lag of 
10** yr (corresponding to an upper Umit to quasar lifetimes) 
implies that s should be larger than ~ 10"^ at the redshift 
z ~ 5. These large tolerances drastically enhance the stochas- 
tic noise. 

- The ratio of the number of topological ghosts to the total 
number of objects (A') is large here, ~ 12 - 17%. In reality, 
however, it can be much smaller than this because of a tiny 
number of ghosts, a large number of objects, or the both, so 
the topological signal will be hidden by the stochastic noise: 
roughly speaking, the former scales as A^^ while the latter 
scales as A^*". 

- Each simulation here took several minutes using an ordi- 
nary personal computer Howeover, calculations for a real- 
istic value of N ~ IQ^ takes an extremely long time that is 
roughly proportional to nCt, - N(N - l)(N - 2)16. 



For example, with the choice of e = lO"'' and A^ ~ 60, the 
stochastic noise has an order of 10^', completely hiding the topo- 
logical signal of 10^. The situation gets even worse for a larger 
A^. 

Thus the method will be impractical in realistic situations 
due to small signal to noise ratio. One way to overcome the 
problem may be using a large n > 4, rather than n = 3. For 
a given pair of n-tuplets to be linked by a holonomy, there 
are n - 2 conditions to be satisfied, thus we have n - 2 fil- 
ters. The possibility P{n), such that a non topological 2n-tuplet 
passes through all the filters by chance, is monotonically de- 
creasing with n. Therefore, the stochastic noise P(n) x nC„, can 
be suppressed if we choose a sufficiently large n. However, as 
mentioned in the previous section, the calculation time roughly 
scales as tjC„ - Nl/{N - ny.nl and calculations for such a large 
n are extremely time-consuming. 

Considering these circumstances, we must conclude that it 
is unrealistic to constrain cosmic topology from the celestial 2D 
distribution of objects. Rather, for the practical purpose, we need 
the 3D catalogs of cosmic objects obtained by spectroscopic ob- 
servations. The new crystallographic method introduced by Fujii 
& Yoshii (2011) can be applied to such catalogs, and provides 
us a knowledge about the shape of the Universe. This will more 
than ever strengthen the motivation for promoting a large-scale 
spectroscopy survey of high redshift objects. 




Fig. 1. The i,-histograms for various types of holonomies: trans- 
lation {top left), half-turn corkscrew motion {top right), quarter- 
turn corkscrew motion {mid left), third-turn corkscrew motion 
{mid right), sixth-turn corkscrew motion {lower left), and glide 
reflection {lower right). Bumps constituted by topological copies 
for these types are seen in the large s, regions for multiconnected 
spaces. 
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Holonomy type 


Number of 
quasars 


Number of 
ghosts 


Topological 
index 5m..if 

u mull 


Topological 
index S cimn 


Translation 


57 


14 


224 


116 


Half-turn corkscrew motion 


60 


20 


733 


175 


Quarter-turn corkscrew motion 


58 


16 


186 


138 


Third-turn corkscrew motion 


57 


14 


318 


235 


Sixth-turn corkscrew motion 


59 


18 


372 


291 


Glide rcllcclioii 


58 


16 


359 


153 



Table 1. Results for various types of holonomies. Topological index S muit is the number of sextuplets satisfying \AD - BC\ < 10 ' 
for a multiconnected space, and S simp is for the simply connected one. 
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